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SUMMARY 


A method has been developed for converting a matrix of direction cosines 
into an equivalent Euler sequence. For a desired rotational sequence, the 
analysis produced a set of five equations which require a sequential calcula- 
tion of the Euler angles. The first angle is used to compute the second and 
the second is used to compute the third. For a given direction cosine matrix, 
the equations show that there are two Euler angle sets which will generate 
that matrix and can he considered to he the result of either a "positive" or a 
"negative" initial rotation. The results of the analysis will work throughout 
the entire 360° angle range and also for the singular cases. For singular 
cases the equations become indeterminate hut the problem can be resolved if 
one recognizes that the two angles not involved in the singular condition 
simply add directly and that only their sum affects the direction cosines. 

The results of the analysis have been generalized into two sets of equa- 
tions, one of which applies to the classical or repeating sequences and the 
other to the nonclassical or nonrepeating sequences. These equations have 
been written as a Fortran TV subroutine and are presented in the appendixes of 
the report . 


INTRODUCTION 


In a research program for a satellite attitude control system it was 
found advantageous to write the equations of motion with the direction cosine 
matrix as the kinematic variable. If parameters such as Euler angles are used 
for these variables, as is often done, it is necessary to be able to convert 
the direction cosine output to Euler angles and the reverse. The conversion 
of Euler angles to direction cosines is simple and can be performed by the 
multiplication of elementary rotation matrices or by the use of a standard 
(e.g., refs. 1 , 2 , and 3 ) which gives the direction cosine matrix in terms of 
Euler angles. It is clear that the calculation of the Euler angles from the 
direction cosines is also possible, but a general method was by no means 
obvious from inspection of the equations involved and no reference could be 
found that gave a general method. In this report, a technique is given for 
performing the conversion. The results are given such that if a direction 
cosine matrix is specified and a rotation sequence given, the Euler angle 


sequence that will produce the direction cosines can be computed. The results 
for the technique have been written as a Fortran IV program presented in the 
appendix . 


TABLE OF SYMBOLS 


A 

E3 (0 1) 

5 2 

0 


Analysis Section 

element of the direction cosine matrix 
direction cosine matrix, 3 x 3 

typical elementary rotation matrix; first rotation angle 
01 , about axis 3 

eigenvector of the rotation subscripted 
rotation angle 


Subscripts and Superscripts 

i, 2>3 coordinate axis for the rotation when used as a subscript 

for E; order of rotation of the angles when used as a 
subscript for 0 

t transpose of matrix 

I,J,K rotation sequence 

th( 1 ) 

I,J,K 
A(I,J) 

EULANG 

ARTN(X(I),Y(I)) 


Computer Program 

angle 0 j_ 
rotation sequence 

a ij 

conversion subroutine name 
arctangent function routine, tan - 


-Vx(i) 
.Y(I) 


ANALYSIS 


The equations for converting the direction cosines to Euler angles will 
be derived for a specific Euler sequence and will then be generalized for any 
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sequence. For a l,ci ,3 rotation sequence, the direction cosines matrix can be 
generated by the multiplication of three elementary rotation matrices. 


A = E 3 ( 0 3 )E 2 ( 0 2 )Ei( 0 i) 


The elementary rotations are given by 

1 0 

Ei(0i) = [ 0 cos 0i 

,0 -sin 0i 



cos 0 2 0 -sin 0 



e 2 (0 2 ) = y 0 1 

sin 0 2 0 


cos d 3 sin 0 3 0 

E 3 ( 0 3 ) = [ -sin 0 3 cos 0 3 0 

0 0 1 


Equation (l) may be written in the following form 

A E^fei) = E 2 ( 0 2 ) 


J 


( 1 ) 


( 2 ) 


( 3 ) 


The coordinates of the eigenvector of the middle rotation are given by 


6 2 - 



w 


If both sides of equation (3) are multiplied by 5 2 , the following equation 
results 


E 3 t ( 0 3 ) A Ei 11 (61)82 = E 2 ( 0 2 ) 6 2 = 6 2 


( 5 ) 


Equation (5) may "be rewritten in the form 


A E i^(0 1)5 2 “ 23(63)62 (6) 

If the following two definitions are made for the portions of equation (6), 

x(0i) = E i t (e 1 )6 2 = cos 0! ) ( 7 ) 

\s in ©i 


sin 03' 

Y(0 3 ) = E 3 (0 3 )& 2 = 1 cos 0 3 

0 


(8) 


Then, equation (6) may be written in the form below, where 0 2 has been 
eliminated from consideration. 


A X(0i) = Y(0 3 ) 


(9) 


Equation (9) is a set of three equations 


a 32 cos 0i+ a 33 sin 0 1 = 0 


( 10 ) 


a 22 cos 0i + ai3 sin 0i = sin 0 3 
a ?? cos 0i + a2 3 sin ©i = cos 0 3 


(ID 


And from equation (3)> by expanding the left-hand side and equating the an 
and a i3 terms to the corresponding terms on the right-hand side, we obtain 


an cos 0 3 - a 2 i sin 0 3 = cos ©A 


a 3 i = sin 0 2 


( 12 ) 


The desired Euler angles can be calculated from equations (10) through (12) . 

It is necessary first to solve for 0i from equation (10), then for 0 3 from 
equation (ll) , and finally for 0 2 from equation (12) . Note that all nine 
elements are used . 


The solution of equations (10) through (12) results in two sets of Euler 
angles because the solution of equation (10) is double valued. For each value 
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of 9 i, unique values of 0 3 and 0 2 are found from equations (ll) and (12) so 
that two distinct Euler sequences result for a given direction cosine matrix. 
Physically the result means that it is possible to reach a given position by 
• starting with either a "positive" or a "negative" rotation 0i followed by 
the appropriate rotations, 0 2 a nd 03 > to give the final position. For a non- 
classical sequence such as the one used in the previous analysis (l,2,3), the 
' relation between the two sequences is as shown below. 


case 1 

0i , 

0 2 , 

^3 

case 2 

01 ± « , 

— 0 2 i jt , 

+1 

CO 

CD 


(13) 


For the classical or repeating type sequence, the relation between the angles 
is different and can be shown to be of the following form: 


case 1 

0i > 

CM 

CD 

03 

case 2 

0i ± it , 

- 0 2 , 

03 


( 14 ) 


For each Euler sequence there is a singular case which requires special 
handling. The repeating sequence is singular when sine 0 2 is zero and the 
nonrepeating sequence is singular when cosine 0 2 is zero. When the 1,2,3 
sequence is singular, the direction cosines a 32 and a 33 are both zero and 
the computation of 9 i is indeterminate. The problem can be resolved by 
noting that 9 i and 0 3 add directly. Since the direction cosine matrix only 
specifies the final position, 0i and 0 3 may have any value so long as their 
sum is correct. In the computation it is necessary to set 0i to some value, 
possibly zero, and then proceed. 


Once equations (10) through (12) have been derived, it is easy to see how 
they could have been obtained by direct inspection of the direction cosine 


matrix 

matrix 

written 

is. 

in terms 

an 

a ip 

a 13 


a 2 i 

a 22 

a 23 


a 31 

a 32 

a 33 


1 

cos 0 2 

cos 0 

3 


-cos 0 2 sin 0 3 
sin 0 2 


cos 0 i sin 0 3 
+ sin 0i sin 0 2 cos 0 3 

cos ©i cos 0 3 
- sin 0i sin 0 2 sin 0 3 

-sin 9 i cos 0 2 


sin 0i sin 0 3 
- cos 0i sin 0 2 cos 0 3 

sin 0i cos 0 3 
+ cos 0i sin 0 2 sin 0 3 

cos 0i cos 0 2 


(18) 
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The following operations 

a 32 = -sin 0i cos 0 2 , multiply by cos 0i 

a 33 = cos 0i cos 0 2 , multiply by sin 0i 

and adding the two equations yields equation (10) 

a 32 cos 0i+ a 33 sin 0 i = 0 


Next, 


a 3 i = sin 0 2 


together with the following operations 

a i 3 = cos 0 2 cos 0 3 , multiply by cos 0 3 

a 2 i = -cos 0 2 sin 0 3 , multiply by -sin 0 3 


and adding the last two equations, yields equations (12) . 

an cos 0 3 - a 2 i sin 0 3 = cos 0 2 
a 3 i = sin 0 2 


Finally, to obtain equations (ll), perform the following operations: 

a i 3 = cos 0i sin 0 3 + sin 0i sin 0 2 cos 0 3 , multiply by cos 0i 

ai 3 = sin 0i sin 0 3 - cos 0i sin 0 2 cos 0 3 , multiply by sin 0i 

and the two equations add to get 

an cos 0i + ai 3 sin 0i = sin 0 3 


Similarly, 

a 22 = cos 0i cos 0 3 

a 23 = sin 0i cos 0 3 


sin 0i sin 0 2 sin 0 3 , 
cos 0i sin 0 2 sin 0 3 , 


multiply by cos 0i 
multiply by sin 0i 


yields 


a 22 cos 0i + a 23 sin 0i = cos 0 3 


The resulting equations (24) and (26) are equations (ll) • 


(19)' 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 


(26) 


6 


An inspection of a direction cosine matrix written in terms of any other 
Euler sequence will show that the same procedure can be followed with the 
.proper choice of terms. Also it should be noted that the solution sequence 
could be in reverse order, 0 3 , Q 1} 0 2 instead of 9 9 3 , Q 2 . 

The derivation that resulted in equations (10) through (12) can be gen- 
eralized to apply to any Euler sequence since Ej(0 2 )6j = If the rotation 

sequence is I, J, K, the direction cosine matrix wi ll be given by 


A - Ej£(0 3 )Ej(0 2 )Ej(0 1) (27) 

The first step is to form the eigenvector of the middle rotation, &j, and then 
compute the column vectors 


X(0i) = E I t (0 1 )5 J ' 

Y(0 3 ) = E k (0 3 )6j , 

The equation to be used for the solution of 0j. and 0 3 is then 

A X(0i) = Y(0 3 ) 


(28) 


(29) 


and 0 2 can be computed by use of equation (3) rewritten in the form 

Ej( 0 2 ) = E K t(0 3 ) A E^©!) 


(30) 


A set of five equations similar to equations (10) through (12) will result 
from the expansion of equations (29) and ( 30 ) • The 12 possible sets of equa- 
tions are given in appendix A. 

The 12 sets of equations that result from equations (29) and (30) can be 
reduced to two sets of equations with an appropriate method of indexing the 
direction cosine elements. One of the two sets applies to the classical or 
repeating sequences and the other to the nonclassical or nonrepeating 
sequences. These equations have been written as a Fortran IV subroutine pre- 
sented in appendix B. For computing, the necessary input to the subroutine is 
the array of nine direction cosines and the desired rotational sequence. The 
output will be the three Euler angles . Since there are two Euler sequences 
for a given direction cosine matrix, the program has been set up to compute 
only one. If the other is desired, equations (13) or (l4) can be used for 
the conversion. Also, in the singular case computations, 0i has been set to 
zero . For a sequence of computations this may cause a discontinuity in the 
output. Smoothing may be accomplished by extrapolation of previous values of 
0i and using this as the output value of 0! and subtracting the same amount 
from 0 3 . 
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CONCLUDING REMARKS 


In the preceding analysis a method was developed for converting a matrix • 
of direction cosines into an Euler angle sequence. For a given rotational 
sequence a set of five equations can he written which call for a sequential 
calculation of the three Euler angles. The first angle is used in the calcu- * 
lation of the second and the second in the calculation of the third. The 
equations show that for a given direction cosine matrix and a specified rota- 
tional sequence there are two Euler angle sequences that satisfy the equa- 
tions. The reason is that the particular values of the direction cosines may 
he produced hy an initially positive rotation followed hy two others of 
correct magnitude or an initially negative rotation and then two additional 
rotations. For the singular case the conversion equations become indetermi- 
nate hut the first and last rotations add directly. The indeterminate condi- 
tions arise because any two values of these angles will satisfy the direction 
cosine matrix as long as their sum is correct. It is necessary to assume some 
value for the first angle (zero or some extrapolated value from previous 
calculations) and proceed with the computations. 

The equations for the conversion have been generalized so that for a 
given direction cosine matrix and any desired rotational sequence the Euler 
angles can he computed. This set of equations has been written as a Fortran 
Fortran IV computer subroutine and is presented in appendix B. 
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APPENDIX A 


EQUATIONS REQUIRED TO CONVERT DIRECTION COSINES TO EULER ANGLES 

FOR ALL 12 EULER SEQUENCES 

1 , 2,3 Sequence 
a33 sin ©i + a 32 cos ©i = 0 
a i3 sin ©i + ai2 cos ©i = sin ©3 
a 23 sin ©1 + a 22 cos ©i = cos ©3 
a3i = sin 0 2 

-a 2i sin 0 3 + an cos 0 3 = cos 0 2 

1 . 3.2 Sequence 
a 22 sin ©i - a 23 cos ©i = 0 

aig sin ©1 - ai3 cos ©i = sin © 3 
-a 32 sin ©i + a 33 cos ©i = cos 0 3 
— a 2 1 = sin © 2 

a 3 i sin © 3 + an cos 0 3 = cos 0 2 

2 . 1.3 Sequence 

-a 33 sin ©i + a 3 i cos ©i = 0 
a 23 sin ©1 - a 2 i cos ©i = sin 0 3 
-a 13 sin ©i + an cos ©i = cos 0 3 
-a 32 = sin 0 2 

a ip sin ©3 + a 22 cos 0 3 = cos © 2 
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2 , 3,1 Sequence 


an sin 0i + ai3 cos 0! = 0 

a 2 i sin 0i + a 23 cos 0i = sin 63 

a 31 sin 0i + a 33 cos 0 X = cos 0 3 

a 22 - si 11 

-a 32 sin 0 3 + a 22 cos 0 3 = cos 0 2 

3 , 1,2 Sequence 
a 22 sin 0i + a 21 cos Q x = 0 

a 32 sin 0i + a 31 cos 0i = sin 0 3 

a^ sin 0i + an cos 0i = cos 0 3 

a 23 = sin 0 2 

-a i3 sin 03 + a 33 cos 0 3 = cos 0 2 

3 . 2.1 Sequence 

an sin 0i - aiz cos 0 X = 0 

a 3 i sin 0i - a 32 cos 0i = sin 0 3 

-a 2i sin 0i + a 22 cos 0i = cos 0 3 

-a 13 = sin 0 2 

a 23 sin 0 3 + a 33 cos 0 3 = cos 0 2 

1 . 2.1 Sequence 

ai3 sin 0i + a^ cos 0i = 0 

a 33 sin 0i + a 32 cos 0i = -sin 0 3 

a 23 sin 0i + a 22 cos 0i = cos 0 3 

a 2i sin 0 3 + a 3 i cos 0 3 = sin 0 2 

a 11 = cos 0 2 
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1,3 >1 Sequence 


a sin 0i - ai 3 cos 0i = 0 

-a 22 sin 0 1 + a 23 cos 0i = sin 0 3 

-a 32 sin 0i + a 33 cos 0i = cos 0 3 

a 3 i sin 0 3 — a 2 i cos 3 i = sin 0 2 

a n = cos 0 2 

2 .1 .2 Sequence 

a 23 sin 0i - a 2 i cos 0i = 0 

-a 33 sin 0i + a 3i cos 0i = sin 0 3 

-a i 3 sin 0i + an cos 0i = cos 0 3 

a i 2 sin 0 3 - a 32 cos 0 3 = sin 0 2 

a 22 = cos 0 2 

2.3.2 Sequence 

a 2 i sin 0i + a 23 cos 0i = 0 

an sin 0i + ai 3 cos 0i = -sin 0 3 

a 3i sin 0i + a 33 cos 0i = cos 0 3 

a 32 sin 0 3 + ai 2 cos 0 3 = sin 0 2 

a 22 = cos 0 2 

3.1.3 Sequence 

a 32 sin 0i + a 3i cos 0i = 0 

a 22 sin 0i + a 2 i cos 0i = -sin 0 3 

a 12 sin 0i + an cos 0i = cos 0 3 

a i 3 sin 0 3 + a 23 cos 0 3 = sin 0 2 

a 33 = cos 0 2 
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3 ,3 Sequence 


a 3 i sin 0i - a 32 
-an sin 0 i + a 22. 
— Si2 1 S in 01+ 3-22 
a 23 sin 0 3 - a 13 
a 33 — cos 02 


cos 

01 - 

0 


cos 

01 - 

sin 0 3 

cos 

01 = 

cos 

03 

COS 0 3 = 

sin 

02 
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APPENDIX B 


FORTRAN IV PROGRAM FOR DETERMINING EULER ANGLES FROM 
A GIVEN DIRECTION COSINE MATRIX 


SUBROUTINE EULANG (l,J,K,A,TH) 
DIMENSION X(3) , Y(3), TH(3), A(3,3) 
IF (I.EQ.K) GO TO 103 
L = I - M0D(J,3) 

IF (L.EQ.2) L = -1 
C = L 

X(l) = A(K,j)*C 

Y(l) = A(K,K) 

TH(l) = ARTN(x(l),Y(l)) 

X(l) = sin(th(i) ) 

Y(l) = COS(TH(l) ) 

101 X(3) = A(l,K)*X(l) - A(l,j)*C*Y(l) 

Y(3) = A(J,J)*Y(1) - A(«J,K)*C*X(l) 

X(2) = -A(K,l)*C 

Y(2) = A(I,I)*Y(3) + A(J,I)*C*X(3) 
GO TO 104 

103 N = 6 - (K + J) 

X(l) = A(K,J) 

L = N - M0D(I,3) 

IF (L.EQ.2) L = -1 
C = L 

Y(l) = A(K,N)*C 
TH(l) = ARTN(x(l),Y(l)) 
x(l) = sin(th(i) ) 

Y(l) = cos(th(i) ) 

102 x(3) = -A(N,N)*X(l) + A(N,J)*C*Y(l) 
Y(3) = a(j,j)*y(i) - A(J,N)*C*X(1) 
X(2) = A(J,I)*X(3) - A(N,l)*C*Y(3) 
Y(2) = A(K,K) 

104 th( 3 ) = artn(x( 3 ),y( 3 )) 

TH(2) = ARTN(X(2),Y(2)) 

RETURN 

END 
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